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1 Introduction 

This is a report on numerical work on extremal Kahler metrics on toric surfaces, using the ideas developed 
in [3pl . These are special Riemannian metrics on manifolds of 4 real, or 2 complex dimensions. They can be 
regarded as solutions of a fourth order partial differential equation for a convex function on a polygon. We 
refer to the preceding article [4j for further background. Motivation for this kind of work can be found in at 
least three different sources: 

• Great effort is expended on the proof of abstract existence theorems for such metrics. Only very rarely 
does one expect to find explicit formulae for the solutions, ft is interesting to explore the geometry of 
the solutions using numerical methods. 

• ft is conceivable that specific numerical results of this general nature could be relevant to physics [7] . 

• The numerical work illuminates the analytical techniques and difficulties involved in the abstract 
existence theory. This is true both of the approximate numerical solutions and of the approximation 
scheme itself, which is founded on a link between algebraic geometry and asymptotic analysis. 

The results described here are obtained using a computer programme which can be applied in principle 
to any polygon and in practice works effectively for a reasonable range of cases. We focus on four such 
cases: particular p-sided polygons with p = 5, 6, 7, 8, which correspond to complex surfaces constructed by 
"blowing up" the projective plane p — 2 times. The status of the rigorous existence theory is different in the 
four cases: 

• For the hexagon, the metric we seek is a Kahler-Einstein metric whose existence was proved by Tian, 
Yau and Siu in the 1980's. Using a rather different approach, Doran et al [6] also recently found 
numerical approximations to this metric. 

• For the pentagon the metric we seek was proved to exist very recently by Chen, LeBrun and Weber [TJ 

• For the octagon, the metric we seek has constant scalar curvature and the abstract existence follows 
from recent work of the second author. 

• For the heptagon, the metric we seek is extremal but not constant scalar curvature, and as far as we 
know is not at present covered by any general existence theorem. 

But we emphasise that these cases are just a small sample of the ones that can be treated successfully, and 
there are many other cases of the fourth kind, where the numerical results outrun, at present, the abstract 
existence theory. 

In each of these four cases we give a selection of geometrical results: again these form a small sample of 
what one can do. The pentagon is particularly interesting because, as Chen, Le Brun and Weber showed, 
it yields a new (conformally-Kahler) solution to the Riemannian Einstein equations. (More precisely, the 



1 Although the simpler results described here were obtained by the second author in 2004, before those of [3]. 
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extremal metric we approximate is very close to another extremal metric which is conformal to an Einstein 
metric.) We verify this numerically, from which one could go on to obtain information about the Einstein 
metric. 

In the final section, we make some comments about ways in which these techniques might be taken further 
in the future. 

We are grateful to Julien Keller for many discussions about this work, and much crucial technical assis- 
tance. In particular, the C++ version of the programme to find geodesies was written by Keller. 

2 The set-up 

2.1 Algebraic metrics 

Our starting point is an integral Delzant polygon P, contained in the square of side length k. So for example, 
this is a picture of our hexagon, with k — 6. 

,(6,6) 



-1 



-1 

-1 

(0,0) (3,0) 

When we write, for example, "the hexagon with k = 12" , we mean the same figure scaled up by a factor 
of 2. So, for the hexagon, we can take k to be any even number (because we need to work with polygons 
having integral vertices). This use of language differs slightly from that in [4], but we hope the meaning 
will be clear. We also indicate in the diagram the self-intersection numbers of the 2-spheres in the manifold 
corresponding to the edges; in this case all — 1. 

Next we consider a collection of positive numbers (a„) indexed by the lattice points v in the closed 
polygon P. Then we define a "Kahler potential" 

= log (l>^- 

This a convex function (of a variable t € R 2 ) whose derivative maps onto the polygon P. The Kahler 
potential defines an "algebraic metric" on the toric surface X corresponding to P. The Lcgendre transform 
of (f> is a convex function u on P. The scalar curvature S is given by the expression S(u) = — ■ q™q x ■ If 
u were to satisfy the partial differential equation 

S{u) = -A, (1) 

for an affine linear function A, then our metric would be extremal. This occurs in two very special cases: 
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• If P is the mi x 777,2 rectangle, with origin as one vertex, then we can take the product of binomial 
coefficients 

— ( mi \ ( 1712 

for integers < v\ < m\, < v 2 < m 2 - 

• If P is the triangle with vertices (0, 0), (fc, 0), (0, fc) we can take the multinomial co-efficents 

fc! 



a 



v 1,1*2 



vi\v-i\{k — v\ — v 2 )V 



The metrics we describe are, respectively, the product of round metrics on S 2 x S 2 and the Fubini-Study 
metric on CP 2 . These are extremal metrics but of course we have not really gained anything new. In general, 
we cannot hope to find an extremal metric exactly, in this way, but we can hope to find good approximations, 
and we can hope that the approximations improve when we increase fc. For example, this is a collection 
of co-efficients a v which give a good approximation to the extremal metric in the case of the hexagon, for 
k = 6: 



2.06 7.17 7.17 2.06 



7.17 38.9 63.5 38.9 7.17 



7.17 63.5 163.6 163.6 63.5 7.17 



2.06 38.9 163.6 257.2 163.6 38.9 2.06 



7.17 63.5 163.6 163.6 63.5 7.1; 



7.17 38.9 63.5 38.9 7.1/ 



2.06 7.17 7.17 2.06 



In this case there is a symmetry under a dihedral group of order 12, so we can reduce to only 6 independent 
coefficients. For clarity we have just written these to a few significant figures. When fc is larger the data 
(a v ) is best kept confined to a computer file, and we will just record the maximum and minimum of the 
a„. It may be helpful to recall here from [5] that the numbers (— loga„) can be thought of, roughly, as a 
"discrete approximation" to the symplectic potential u. Saying much the same thing in a different way, if 
we view u as a function on P depending on parameters a v , then at a given point x G P the value of u is 
strongly influenced by the a„ for v close to x, but only little by those for far away v. These somewhat vague 
statements become more precisely true when the parameter fc is large. 

If a is any affine-linear function then changing the coefficients a v to e a ^a v does not change the metric 
described by the data. We can normalise the a v to take account of this, as follows. Let {p\,P2) be the centre 
of mass of the polygon. We arrange that 

log a„ = , ^(i>i - pi) loga„ = , ^(^2 - P2) log a v = 0. (2) 
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2.2 Decomposition of the curvature tensor 

Here we want to extend slightly the discussion of the curvature tensor of a toric manifold in 0], for the 
two-dimensional case. 

We begin in a more general setting and consider a complex Kahler surface Z with an anti-holomorphic 
isometric involution a fixing a submanifold E of two real dimensions. Thus S is a "real form" of Z . We 
want to understand the curvature tensor of Z at a point p of E. We have TZ p — TE p N p where N is 
the normal bundle. Fix an orientation of E, at least locally. This determines an orientation of the normal 
bundle so, using the metric, we can define complex structures on each of TS and N, making them complex 
line bundles. With suitable choice of signs the resulting complex structure J_ on TX P induces the opposite 
orientation on Z from that of the Kahler structure. Thus we have a unit anti-self-dual 2-form at p and 
we can write 

A p = R6»_ L, 

say, where L is a two-dimensional oriented real vector space. Thinking of TZ p as the complexification of 
TS p we see that there is a natural map from the trace-free symmetric 2-tensors Sq = SqT* to the trace-free 
Hermitian forms Aj' 1 on Z at p. In a standard basis this is just the obvious map which takes a symmetric 
real matrix to a Hermitian matrix. It is easy to check that the image of this is L, so we can identify L with 

s 2 

Now we use the well-known decomposition of the curvature tensors of Riemannian 4-manifolds, and in 
particular Kahler surfaces. In the general Riemannian case, the curvature tensor, viewed as a symmetric 
element of A 2 A 2 , decomposes into four pieces according to the ±-self-dual splitting of the 2-forms, and 
two of these are equivalent by the symmetry. The A + ® A~ block is equivalent to the trace-free Ricci tensor. 
When the 4-manifold is Kahler the A + A + piece (the self-dual Weyl tensor W + ) is entirely determined by 
the scalar curvature. The scalar curvature also determines the trace of the piece W~ in A~ (gi A - . Further, 
since the Ricci tensor has type (1,1), the A + (g) A - piece lies in (uj) (g) A~ , where uj is the Kahler form (which 
is self-dual). 

Putting these observations together, we see that the curvature tensor of Z at p has at most the following 
six irreducible components 

1. The scalar curvature S. 

2. A component of the trace-free Ricci tensor in (uj) ®Lc A + ® A - . This can be regarded as an element 
p of s 2 . 

3. A component of the trace-free Ricci tensor in the 1-dimensional space (uj) ® (6L) C A + <g> A - . 

4. A component of the anti-self-dual Weyl tensor W- in ® (#_) C A~ <E> A - . 

5. A component of W~ in (0_) ®LcA+® A+. 

6. A component w of W~ in the 2-dimensional space Sq(s 2 ). 

It is an exercise to check that, since the curvature tensor is preserved by the action of a on TZ p , which 
reverse the Kahler form, the components (3) and (5) vanish. We can identify the component p in (2), which 
determines the trace-free Ricci tensor, with a quadratic differential on E — an element of T*S <8>c T*S. 
Likewise, we can identify the component w in (6) with a quartic differential on E. The scalar component 
(4) can be written as S — 6K where K is the Gauss curvature of the induced metric on E. The conclusion 
is that the curvature tensor is determined by S, p, K, w. 

We now apply this to our toric surface X, taking E to be the real form Xr. We work in the local co- 
ordinates £i, £2) but the point of the above discussion is to see that the tensors we obtain extend smoothly over 
this compact surface. (If we want to use the complex description, with quadratic and quartic differentials, we 
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should go to the oriented cover.) The four-dimensional Riemannian metric associated to the convex function 
4>{h,t2) is given by 

<pabdt a dt b + cf) ab d9ad9 b . 

Here 9±, 9i are angular coordinates and the 4> ab are the second derivatives of <p with respect to the t- variables. 
The 9 a coordinates never appear in the calculations and, because of the torus action, the curvature tensor 
over the whole of X is obviously determined by its restriction to Xr,. 
We define a tensor 

Fgbcd — <Pabcd - ^ 4>abi(f>cd,j<p lj , (3) 
ij 

where (</> y ) is the inverse of the Hessian V 2 </> = (<fi a b)- The tensor F is symmetric in the pairs of indices (ab) 
and (cd) and also under interchange of the two pairs. The scalar curvature is 

s = ]T F abcd r b t cd . 

abed 

The quadratic differential p corresponds to the trace-free Ricci tensor and is given by 

Pab = F abcd (j) Cd — —(^ab- 

ed 

The other scalar component is determined by the Gauss curvature K of the metric restricted to the (t\, t?) 
variables, which is 

K = (F1122- F 1212 ) det{(j) ab ) 

The final component w, in the space of quartic differentials, is obtained by projecting F to the the kernel 
of these contractions. We will call this the "Weyl component" , although strictly the Weyl curvature of the 
4-manifold also contains contributions from S and K. With our conventions we have 

| Riem | 2 = is 2 + 1(S - %Kf + | W | 2 + \\P\ 2 - (4) 

Examples 

• If X is CP 2 , with the standard Fubini-Study metric scaled to S = 1, then 

5 = 1 p = \w\ 2 =0 K =1/6 (5) 

and |Riem| 2 = 1/3. 

• If X is 5 2 x 5 2 , with each factor having the round metric of the same area then, scaled to S = 1, 

5=1 p = M 2 = l/8 K = (6) 
and |Riem| 2 = 1/3 + 1/24+1/8 = 1/2. 
We note also that Chern-Weil theory gives a global relation 

/ |Riem| 2 -5 2 = (2^) 2 (p-6), (7) 
Jx 

where p is the number of vertices of the polygon. 

We will also consider the "Bach tensor" . According to Chen, LeBrun and Weber, p], this can be written, 
on any Kahler surface, as the trace-free part of 

5 Ricci + iddS, 
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where S is the scalar curvature and Ricci is the Ricci form. If the metric is extremal this is a harmonic (1,1) 
form. In our setting we can write the Bach tensor as a quadratic differential given by 

Bab = Spab + S ab ~ (<% ^ )0ab , 

where S a b denotes the partial derivative g f g t . If the scalar curvature S is positive and we consider the 
conformally related metric g = S~ 1 g then the trace-free Ricci tensor of g is SB a i,. In particular if the Bach 
tensor vanishes then this conformal metric is Einstein. All this is explained in [1]. 

Suppose we have any trace- free symmetric 2-tensor 9 ao on Xr. This defines a (1,1) form 9 a bdT a dr b on 
X. The condition that this form is harmonic is simply that 

dOgb _ dO ac 
dt c dt b ' 

This equation is similar to that defining the holomorphic quadratic differentials on Xr, in terms of the 
conformal structure on Xr. More precisely, the equation can be written as 89 + L{9) — 0, where L is 
a certain algebraic operator, determined by <f>. Let us call the solutions of this equation "0-holomorphic 
quadratic differentials" . Then we see that for an extremal metric the Bach tensor can be viewed as a <fi- 
holomorphic quadratic differential on Xr. One can show further that, if the Bach tensor vanishes, then 
the product Sw is a "0-holomorphic quartic differential" on Xr, in a similar sense. Just as for genuine 
holomorphic differentials, the zeros of </>-holomorphic differentials occur with positive signs, so the number 
of zeros, counted with multiplicity, is fixed by topology. 

2.3 Integration 

The foundation of our approximation scheme is the numerical integration of functions on the manifold X. 
Since the functions will all be invariant under the torus action the integration is really in 2 dimensions 
rather than 4, but the presence of the fixed points of the action, corresponding to the boundary of the 
polygon, means that this reduction is not completely straightforward. Suppose given a coefficent set (a u ). 
The functions we want to integrate are given explicitly, by standard formulae, in the t coordinates, so one 
approach would be to fix a lattice covering a large square in the (t%, i2)-plane and approximate the integrals 
by the corresponding sums. There are two drawbacks to this. First, we will get errors from the contributions 
to the integrals outside the square. Second, in these t co-ordinates the calculations will involve ratios with 
very small numerator and denominator when |t| is large, which leads to numerical errors. Of course this is 
a reflection of the fact that these co-ordinates are not valid over the whole of X. 

To get around these difficulties we proceed as follows. Suppose we have a fixed convex function U$ on 
the polygon P, computable in the ^-coordinates. The derivative of Uq gives a map DUo : P — > R 2 . On the 
other hand the derivative of the Kahler potential <j>, defined by the (a„) gives a map D<f> : R 2 — > P. So we 
have a composite map Dcf> o DUq : P — > P. Let J(x) be the determinant of the derivative of this map. Fix 
a small lattice spacing h = k/N and let xj a ^ run over the points in the intersection P n N~ 1 Z 2 . Let / be a 
T 2 -invariant function on X, regarded as a function fit) on R 2 . The integral of / over X is approximated 
by the sum 

h 2 J2 J fe M ) f(DU (x^)). (8) 

a 

If we think of / as a function on P, then we are approximating the integral by a weighted average of the 
values at the array of points x^ — Dcf> o DUq(x^). Recall that we have a system of charts covering X, 
each chart associated to a vertex of P. For each point a;^ Q ^ we choose a suitable chart containing the point 
DUq(x^). Thus, after a linear transformation, we can suppose that this vertex is the origin and that P 
coincides, near the origin, with the standard model {xi > 0}. The corresponding chart is formed by taking 
complex co-ordinates 

z a = exp(-(t„ + i9 a )), 
and the chart is chosen so that we work at a point with \z a \ < 1. In these coordinates 



G 



4> = Y, a *M 2vi \ z *\ 2v '- 

If 

We perform all our calculations in these coordinates to evaluate the contribution to the sum (8) from the 
point xj a \ 

We choose the function Uo to be buo where uo is Guillcmin's "admissible symplcctic potential" 

U = J^(A r - Cr) log(A r - C r ), 
r 

in the notation of [4], and b is a positive constant. If we took b = 1 the pattern of points x^ would be similar 
to the original pattern x^ . Indeed, if it happened that <f> were the Legendre transform of uq then these sets 
would be identical. With this choice, the integration scheme would have errors of order iV _1 arising from 
the boundary of P. So we take b > 1, which has the effect of increasing the density of the set x^ near the 
boundary and reducing the error term. In fact we have taken 6 = 2 which seems to give sufficient accuracy. 
(We test the scheme by calculating integrals whose value we know exactly, such as the total volume of the 
manifold or Chern-Weil integrals: for a reasonable number of integration points we typically get errors in 
the range 1CP 6 — 1CP 4 .) The whole procedure is illustrated in Figure 1. The hexagon on the left shows the 
lattice of points 2 in P. Those that map into a particular chart, corresponding to the top left vertex, 
appear in the shaded region in the left hand diagram. This chart is represented by the small diagram in the 
middle, with co-ordinates in the square corresponding to \z\\, \z%\. The points where we calculate, in this 
chart, appear as dots in the square. The hexagon on the right displays the set of points x^ a \ with those 
corresponding to the given chart shaded. 



Figure 1: The integration scheme 




3 Numerical algorithms: balanced metrics and refined approxi- 
mations 

Our problem is to find coefficient sets {a v ) such that the corresponding metric is close to an extremal metric. 
We consider two different conditions on the (a v ) which we can hope have exact solutions, close to the extremal 
metric. The first is the "balanced" condition. The general context for this involves the L 2 norm on the spaces 
of holomorphic sections of line bundles, but in this toric situation we can write everything out explicitly, 
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ignoring the algebra-geometric interpretation if we please. Given (a„) we write 



D(t) = 



E 



do cf>(t) — log D(t). For each lattice point /iePwe have a function 



D(t) 



Then we perform the integrals 



h = [ U(t) detV 2 4- 

JR 2 



We have written these as integrals over R 2 but of course we evaluate them numerically using the scheme 
described above. Now set 

at = I-\ (9) 

We say that the coefficient set (a„), or the metric it defines, is balanced if (a*) is equivalent to (a„) in the 
sense that, for some affine linear function a, 

a* = e a ^a v . 

General theorjH tells us that the solutions of this balanced condition are essentially unique, and that if an 
extremal metric exists then balanced metrics exist for large enough k and converge to the extremal metric 
as k — > oo. Moreover there is a simple algorithm for finding these balanced metrics. Recall that we fixed 
a way (2) to normalise the coefficients to take account of the action of the affine linear functions. Let Ti. 
denote the set of normalised coefficients. We take any set of normalised co-efficients (a„) and compute a*. 
Then we normalise the a* to get a new set (a**). This defines a map T : H — > Ti. and the balanced metrics 
are exactly the fixed points of this map. These can be found by starting with any point a = {a v ) in Ti. and 
taking the limit of iterates T n (a) as n — > oo. 

This procedure provides a robust and straightforward way to find balanced metrics numerically, and these 
give reasonable approximations to the extremal metrics. The drawback is that they do not give very accurate 
approximations, for practical values of k. One can only expect the balanced metric to differ by 0(k~ 1 ) from 
the extremal metric. Thus, as in [3] Sec. 2.2.1, we seek to extend the ideas to "refined approximations". 
Recall that, on the polygon P, the extremal metric is a solution of the equation S — A where A is an affinc- 
linear function of the symplectic coordinates Xi, determined by the elementary geometry of the polygon. The 
extremal metric is a minimum of the modified Mabuchi functional T which is a functional on the space of all 
Kahler metrics. By restricting to the algebraic metrics we get a function F on 7i. We say that a coefficient 
set (a„) is a refined approximation if it is an absolute minimiser of F. 

The condition that (a v ) is a critical point of F can be expressed directly, without explicit reference to the 
Mabuchi functional. We regard the scalar curvature S as a function on R 2 , in the obvious way, and write 
A for the composite A o D<j) so, expressed on R 2 , the extremal equation is S = A. For each /u, we form the 
integrals 

r lfl = [ {S- A) f^t) detV 2 <p dt. (10) 

JR 2 

Then (a„) is a critical point for F if and only if all the r)^ vanish. 

The abstract theory of these refined approximations is less clear-cut than for the balanced metrics. But 
it is reasonable to expect that, assuming an extremal metric exists, 

• When k is sufficiently large there is a unique refined approximation and it is close to the balanced 
metric. 



2 Strictly, some of this theory has not been written down for the case at hand, but to simplify exposition we will ignore this 
here. 
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• As k — *■ oo the refined approximation converges to the extremal metric faster than any power k m . 

We refer to the discussion, for a related problem, in [3J, Sections 4,5 and move on to the main issue, for 
our present purposes, which is finding these refined approximations in practice. (Notice that the "refined 
approximation" is an exact concept, and of course we can only hope to find approximations to it.) 

We will give a very sketchy discussion to motivate the procedure for finding refined approximations which 
we use. Imagine we start with the balanced metric and write So for the difference Sq — A, where Sq is the 
scalar curvature of the balanced metric. We compute the "error coefficients" rj = {r) v ) as above. Suppose we 
have a collection of small numbers (e„), which we will refer to as "correction coefficients", labelled by the 
lattice points in P. Then we can modify the map T to a new map T e by changing the definition (9) to 

at = (l + e v )I- x . 

If e is small we can hope that iteration of T e will converge to a fixed point, close to the balanced metric. This 
fixed point has some scalar curvature S± and we set Si = Si — A. One can argue that Si is approximately 
Sq — € vfv Now let the L 2 -projection of Sq onto the span of the functions /„ be J2 a vfv ^ we choose 
= &v then we can hope that Si is small compared with Sq, i.e that the fixed point of T e is a better 
approximation to the extremal metric than the balanced metric. But by definition, 

where 

M IU ,= ( / M a)/,te)det(V 2 0) dt. 
Jr. 2 

So, if we can invert the matrix M we should define the correction coefficients e in terms of the error coefficients 
V by 

e = M-\ 

and iterate T e to find a fixed point which is a better approximation than the balanced metric. 

Now the matrix M will be very large and, even when practical, inverting it exactly will probably not be 
very wise. As explained in [3], Sect. 4. 2, the matrix M can be expected to be an approximation to the heat 
operator operator e _A / fc on the manifold. So we expect that there are many eigenvalues of M which are 
slightly less than 1, but also some very small eigenvalues. (This is confirmed by numerical tests.) So the 
exact inverse of M will be very large. It is better to invert only on the span of the eigenvalues which are 
close to 1. But on this span the inverse is close to M itself. The upshot is that, rather than inverting the 
matrix, we could just set r\ = ce for some appropriate fixed constant c. (In fact, after experimentation, we 
set c = 0.75.) Then we can still hope that the fixed point of T e is a better approximation than the balanced 
metric. 

With this motivation in place we describe the actual procedure we use. At each stage we have a set of 
"correction co-efficients" e v . At the outset these are set to 0. We iterate the map Te until we are close to a 
fixed point. ( Here "closeness" is judged by the criterion that the ratio a** ja v differs from I by at most some 
small number c' which, after experimentation, we set to .0008.) Then we calculate the "error co-efficients" 
r\ v by performing the integrals (10). We modify the correction coefficients to a new set 

= e v + -75^, 

and proceed to iterate Ti , and so on. A fixed point of the whole procedure should correspond to the refined 
approximation, with r\ v = 0. 

In practice, leaving aside the very sketchy rationale for the procedure, this process does converge, for 
practical purposes, in all the examples we have tested. After a very long while, we find experimentally that 
the procedure typically oscillates between two values, but these differ only in the seventh decimal digit so 
effectively we have reached a fixed point. 
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On the negative side, this procedure is painfully slow when k is large. Just as in [3J, we can hope that the 
technique could be improved with a more detailed and sophisticated analysis. 

We graph the L 2 norm of S — A, plotted on a logarithmic scale, under this procedure; for two of the 
cases discussed in the next Section. There is a short initial period in which the decrease is very rapid. Then 
there is long intermediate period where the decrease is approximately exponential, until the norm is quite 
close to its limiting value. After this the decrease is much slower. (We expect that the first period is that in 
which nonlinear effects are important and that in the second period the iteration is close to the iteration of 
a linear map, which in turn should be close to the linearisation of the "Calabi flow" .) 

l.OE+O-j 
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Figure 2: Evolution of L 2 error: the pentagon with k — 20 (above) and the heptagon with k = 45 (below) 
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4 Numerical results 



Suppose we have an algebraic metric which we hope to be a good approximation to an extremal metric. To 
compare across different values of k we rescale so that the average value of S is 1. Setting S = S — A we 
compute three measures of the size of the "error" : 

• The L 2 error. That is to say, the square root of the average of S 2 . 

• The maximum and minimum values of S. 

• The "maximum and minimum normalised values" of S. To motivate this, consider a situation where, 
on average, S is very small compared with the full curvature tensor. Then an error S might appear 
large compared with the average value of S 1 , while still being small compared with the full curvature, 
and hence it terms of the true geometry of the manifold. More generally, we may have a situation 
where the curvature tensor of the manifold is much larger in some regions than others — in other words 
the natural scales on which to study the manifold arc different in different regions. To take account of 
this we define the maximum and minimum normalised error to be the maximum and minimum values 
of , R r[ , , where |Riem| is the pointwise norm of the full curvature tensor. 

4.1 The hexagon 

Recall from the introduction that this is the polygon which corresponds to the Kahler-Einstein metric on 
the blow-up of the projective plane at three general points. In this special case there is another, simpler, 
approach which one can try — the "Tk algorithm" described in [3], and in a separate programme we have 
verified numerically that this does indeed converge. But here we will just treat the manifold by our general 
procedure — designed for extremal metrics. See also the discussion in [8]. 

We find that we get a useful approximation to the Kahler-Einstein metric with the low value k = 6. The 
co-efficients of the refined approximation appear in the figure in Section 2.1, in the table below we record 
the largest and smallest of these ccfficicnts and the "error measures" . 



Max. Cocff. 


257.2 


Min. Cocff. 


2.062 


L 2 error 


.0077 


Error maximum 


.021 


Error minimum 


-.14 


Normalised error maximum 


.013 


Normalised error minimum 


-.058 



Taking larger values of k we can do much better. There seems to be no practical gain in going beyond 
k = 20 which yields: 



Max. a v 


5.52 x 10 4 


Min. a v 


1.58 x 10~ 4 


L 2 error 


1.33 x 10-" 


Error maximum 


1.23 x 10~ 5 


Error minimum 


-5.01 x 10" b 


Normalised error maximum 


5.8 x 10~ B 


Normalised error minimum 


-2.0 x 10" B 
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In the next table we record data about the size of the components of the curvature tensor. 



Max. Riem| 


2.48 


Min. |Riem| 


0.61 


Max K 





Min K 


-0.817 


Max. \w\ 


2.09 


Max. \p\ 


2.7 x lQ- b 



We see that the trace-free Ricci tensor p is indeed very small and effectively the only non-trivial compo- 
nents of the curvature tensor are the "Gauss" and "Weyl" components. The Gauss curvature, times —1/4, is 
plotted in Figure 3. The Gauss curvature vanishes at the centre of the hexagon and is everywhere negative. 
(We know a priori that it has to be negative on average, since the Euler characteristic of 1r is —2.) The 
norm of the Weyl component is plotted in Figure 4. We see that the Weyl component also vanishes at the 
centre of the hexagon. Recall that this Weyl component is a "0-holomorphic" quartic differential on the 
oriented double cover of Xr. The number of zeros, counted with multiplicity on this double cover is 16. 
Since there are 8 copies of P making up this double cover, we can see from symmetry arguments that the 
only possibility is a double zero at the centre of P. It is interesting to compare the curvature tensor of this 
manifold with those of the standard examples in (5), (6). 

Figure 3: -1/4 times the Gauss component of the Curvature, for the hexagon 




4.1.1 Geodesies 

Using our numerical solutions for the Kahlcr-Einstein metric we can compute geodesies in this Riemannian 
manifold. The inner products with the two Killing fields defining the torus action are conserved quantities 
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Figure 4: The norm of the Weyl component of the curvature, for the hexagon 




which we refer to as "angular momenta" . Fixing the angular momenta, the geodesic equations can be 
regarded as equations for a particle moving around the polygon in a potential. We do not take this any 
further here, and include this discussion only to demonstrate that one can effectively study the Riemannian 
geometry of these solutions. One could take the ideas much further. However the solutions display some 
interesting phenomena. When the angular momenta are reasonably large compared with the horizontal 
component of the velocity the solutions behave effectively as an integrable system, essentially confined to 
tori in the phase space. As the angular momenta arc reduced the dynamics changes its nature and the 
solutions are more complicated. Of course, we might expect the latter from the fact that when the angular 
momenta are zero the equations can be viewed as the geodesic equations on the surface Xr and since this 
has negative curvature the system cannot be integrable. 
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Figure 5: A geodesic in the Kahler-Einstein metric, showing approximately integrable behaviour 
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4.2 The pentagon 

Here we discuss the pentagon, shown with k — 2, which corresponds to the blow-up of the complex projective 
plane at two distinct points, with the Kahler class equal to C\. 



-1 




(0,0) 



-1 



(2,0) 



The metric we seek is an extremal metric but not of constant scalar curvature. The numerical results, 
regarding speed of convergence and the size of the errors, are very similar to those for the hexagon. 
With k = 20 we get 



Max. a v 


6.41 x 10 4 


Min. a v 


3.34 x 10~ Y 


L 2 error 


1.02 x 10~ 6 


Error maximum 


.1.7 x 10~ b 


Error minimum 


-6.4 x 10~ 5 


Normalised error maximum 


7.4 x 10" B 


Normalised error minimum 


-2.5 x 10~ 5 



The curvature data is given by 



Max.|Riem| 


2.57 


Min.|Riem| 


0.56 


Max. K 


0.132 


Min.K 


-.906 


Max.|u;| 


2.19 


Max |p| 


.4378 



The curvature functions arc plotted in Figures 6,7,8. (In Figure 8 one can see the approximate location 
of a single zero of the Weyl component, as predicted by topology.) 



4.2.1 The conformally-Kahler Einstein metric 

As we mentioned above, the existence of this extremal metric was proved recently by Chen, LeBrun and 
Weber pp. In fact they established the existence of a family of extremal Kahler metrics on this manifold, 
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Figure 6: The size of the trace-free Ricci tensor, for the pentagon. 
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Figure 7: —1/4 times the Gauss curvature, for the pentagon 




depending on a real parameter a which determines the Kahler class. The class we are considering has a = 1. 
Chen, LeBrun and Weber show that for a particular value do ~ .986 . . . the Bach tensor vanishes and the 
extremal metric is conformal to an Einstein metric, with conformal factor the scalar curvature. This extremal 
metric is not directly accessible with our method, which is confined to rational values of the parameter with 
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Figure 8: The size of the "Weyl component", for the pentagon 




fairly small denominator. However, since ao is close to 1 we can expect that our extremal metric has small 
Bach tensor, leading to a reasonable approximation to the Einstein metric. We verify this numerically. The 
norm of the Bach tensor is shown in Figure 9. It is less than .005 over most of the manifold, and has 
maximum value about .02. This is comparable with (1 — do), as one would expect. It would be possible to 
compute geometric quantities associated to the Einstein metric to this degree of accuracy. 
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Figure 9: The size of the Bach tensor, for the pentagon 
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4.3 The octagon 



We now move on to the octagon depicted, which corresponds to the projective plane blow up at five points, 
or S 2 x S 3 blown up at four points. The symmetry of the polygon shows immediately that the Futaki invriant 
vanishes, so the metric we seek has constant scalar curvature. The existence of this metric follows from the 
general result in [5]. (An argument of Zhou and Zhu [9] shows that this polygon satisfies the "positivity" 
criterion in the general existence theorem.) 




We find that we need to take k significantly larger than in the previous two cases, to get a reasonable 
approximate solution. Taking k = 18 we obtain 



Max. a„ 


77349.5 


Min. a„ 


.00051733 


L 2 error 


.00522 


Error maximum 


.0334 


Error minimum 


-.0918 


Normalised error maximum 


.011 


Normalised error minimum 


-.0135 



When k = 36 we get 



Max. a v . 


2.16 x 10 s 


Min. a v . 


9.94 x 10- 1U 


L 2 error 


1.07 x 10" 4 


Error maximum 


1.32 x 10" 3 


Error minimum 


-3.48 x 10~ a 


Normalised error maximum 


2.69 x 10~ 4 


Normalised error minimum 


-5.07 x 10~ 4 
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The reason why we need to make k larger than for the previous two polygons is clear when we look at 
the size of the curvature tensor: 



Max. Riem| 


6.86 


Min.|Riem| 


.644 


Max. K 





Min.K 


-2.60 


Max. w 


5.87 


Max \p\ 


1.35 



The maximum curvature in this case is nearly 3 times that for the hexagon and pentagon. Now \f~k~ 
behaves as a length-scale parameter in the asymptototic theory, and curvature scales as length -2 . So we 
would expect that for the octagon we should take k about 3 times as large as for the pentagon or hexagon, 
to obtain comparable accuracy. The figures above are consistent with this. The size of the full curvature 
tensor and of the trace-free Ricci component are shown in Figures 10,11. Notice that the curvature is much 
larger near the boundary of the polygon. Moreover it is larger near edges corresponding to the —2 curves 
than the —1 curves. The Chern-Weil formula (7) shows that as the number of edges of the polygon increases 
the full curvature tensor grow, on average, compared to the scalar curvature. Also, the self-intersection 
number of a boundary sphere can be written as an integral of an expression involving the curvature over the 
corresponding edge, so the curvature must become relatively large there. But the actual concentration of 
the curvature is more dramatic than one might have expected on such general grounds. 

Figure 12 shows the size of the "error" S — A over the polygon. As we would expect the error is much 
larger near the boundary, where the curvature is large (which is reflected by the normalised errors in the 
tables above). Notice the intricate oscillatory pattern in the function, which one can see is related to the 
lattice spacing. 



20 



Figure 12: The size of the error S 




4.4 The heptagon 




The metric we seek is an extremal metric, not of constant scalar curvature. The results follow a similar 
pattern to the previous case, only more so. We take k — 45 (which is about as large as seems reasonable 
with the programme), and we obtain: 



Max. a v 


3.88 x 10 y 


Min. a„ 


2.44 x 1CT 21 


L 2 error 


2.6 x 1(T 3 


Error maximum 


.085 


Error minimum 


-.23 


Normalised error maximum 


.015 


Normalised error minimum 


-.016 



The error is concentrated near the boundary and particularly near the —3 curve, which is where the 
curvature is largest. The modulus .23 of the minimum value looks large, and one might think that the 
approximation is not very satisfactory, but the normalised value is much smaller. Thus we believe that this 
does model an exact solution quite well. The maximum value of the curvature is nearly 6 times that for the 
hexagon, so we should expect to achieve an accuracy similar to that for the hexagon with k <~ 45/6, and the 
results are in line with this. Notice that the curvature function in Figure 13 is almost constant, at a value 
about 0.7, over most of the polygon. 



Max. Riem 


13.9 


Min.|Riem| 


.537 


Max. K 


.249 


Min.K 


-5.48 


Max. \w\ 


11.7 


Max. \p\ 


4.27 
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Figure 13: |Riem| for the heptagon 
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Conclusions 



Rate of convergence 

As we have stated in Section 3, one can expect that the refined approximations converge to the ex- 
tremal metric faster than any inverse power. But, even if proved rigorously, this is only an asymptotic 
statement and does not give any information about particular values of k. To examine this experi- 
mentally, we consider the L 2 error as a measure of the distance to the true solution, working with the 
hexagon, over a range of different values of k: 







k 


4 


6 


8 


10 


12 


16 


20 


Error 


.019 


.0077 


.0012 


.00052 


.00010 


8.0 x 10~ y 


1.3 x 10~ B 



The data fits reasonably well with an exponential decrease 10 ( k /^ 1 . 

Speed and size of the calculation For a fixed "shape" of polygon the number of elementary calculations 
required in each step of the iteration is proportional to N 2 k 2 , since we have to evaluate 0(k 2 ) functions, 
and their derivatives at 0(N 2 ) points. In the more general case of a d-dimensional toric variety this 
would be N d k d . As k increases we need to increase N, since the functions we are integrating become 
increasingly localised. Simple asymptotics suggests that one should take N = 0(Vk), so the number 
of elementary calculations is 0(k 3d ^ 2 ). However, at a given point, most of the functions are extremely 
small, so one can hope to reduce the calculation by ignoring them. This idea was implemented to 
some extent in our programme, which "truncates" the sums to ignore the very small terms. This does 
not affect the accuracy noticeably but the gain we have achieved so far is only modest. With a more 
sophisticated approach, one could hope to make substantial gains, particularly when k is large, or 
in higher dimensional problems. Simple asymptotics suggests that one might be able to reduce the 
calculation at each point to 0(k d ' 2 ), giving 0(k d ) for each iteration step. 

The convergence is also slower as k increases. One expects that the number of steps required to 
effectively reach a fixed point will be 0(k 2 ). (For example, in Figure 2 the number of steps to reach 
the marked "corner" in the graph for the pentagon, with k — 20, is about 500, while for the heptagon, 
with k = 45, it takes about 2200 steps: then 2200/500 — (45/20) 2 .) So in sum we can hope that 
the length of time needed to run the whole procedure, in dimension d, will be 0(k 2+3d ) without the 
truncation of the sums described above, and perhaps 0(k 2+2d ), with this truncation. 

Scope for improvement 

The most striking feature of our results is the extent to which the curvature of the solutions becomes 
concentrated around the boundary of the polygon or, in terms of the 4-manifolds, around the excep- 
tional curves. This is particularly true of the curves with large negative self- intersection. The obvious 
way to handle this efficiently would be to adopt a "multiscale" technique. This idea can be related to 
the "Veronese embedding" in algebraic geometry, as discussed in [4|. Suppose we have an algebraic 
metric cu defined by an array of coefficients a„, for a certain lattice-size ko- The Kahler potential is 
4> = log(^ a v e vt ). Then we can write 

0=ilog(£> e -) 2 )- 

We can expand out the sum to write 
where 

B ^ ^ Q>V\ Ojjy 2 . 



26 



This means that, after rescaling, we can obtain the same metric ui using the array B but with lattice 
scale 2fco. Now suppose we have chosen a subset G c 2k$P n Z 2 . Then we can consider varying our 
array to 

where c M is zero if /i is not in G. The potential function is 

log^^e^ + E^')- 

If the number of elements of G is relatively small then we can calculate this function and its derivatives 
relatively quickly. But if we choose G appropriately we can hope to get a much more accurate approx- 
imation to the true metric over the high curvature regions. For example in the case of the heptagon we 
might take fco = 30, which should suffice to describe the metric over the interior of the polygon, and a 
set G consisting of lattice points near to the boundary, so that we would effectively be approximating 
the metrics in that region as though we were working with k = 60. 

Of course we can imagine iterating this construction to a yet larger lattice scale 4fc , and so on. Also, 
there is nothing special about the factor 2. But what is needed to develop this approach is to find 
some reliable algorithms for updating the different coefficients. 
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